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Abstract: Female inflorescence of Betula platyphylla was sampled at an interval of each two days to analyze the background of gene ex¬ 
pression in floral phase. On the basis of SMART strategy, the driver cDNA was obtained from total RNA of the last sample and the tester 
cDNA was from that of the others by RT PCR which were subsequently used to construct a subtracted cDNA library. The result of the ESTs 
(expression sequence tags) blastX showed that the genes in the subtracted cDNA library could be mainly clustered into 5 groups related to 
metabolism, transportation and signal transduction, cell cycle, stress response, and regulation. The relationship between gene expression and 
development was also discussed. 
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Introduction 

Genetic and molecular studies have uncovered a large number 
of genes that control different steps in flower development in¬ 
cluding flowering time, flower meristem, and organ identity. In 
particular, the ABC model has been proposed for the specifica¬ 
tion of floral organ identity (Coen and Meyerowitz 1991), and 
the discovery of D and E function genes perfected the regulatory 
network about floral development (Colombo et al. 1995; Pelaz et 
al. 2000). At present a great advance has been made to elucidate 
the genetic mechanism of floral development in many species. 

Birch is a kind of eminent tree species in architecture, furni¬ 
ture, and paper making. At present a breakthrough has been made 
in accelerating genetic breeding by the way of manpower early 
flowering in which process the physiological information was 
monitored (Su etal. 2000). However, few of researches about the 
gene expression pattern of flowering were reported, and until 
recently a few genes were found for regulating floral develop¬ 
ment in Betula pendula which have the function as API (Elo et 
al. 2001), PI (Jarvinen et al. 2003), SEP , AC (Lemmetyinen et al. 
2004) in Arabidopsis or SBP (Lannenpaa et al. 2004) in Snap¬ 
dragon. In this study, we report some new study in floral devel¬ 
opment process of Betula platyphylla, which will found a basis 
for future work of cloning genes by the RACE (rapid amplifica¬ 
tion of cDNA ends) method. 

Material and method 

Female inflorescence was sampled from birch forest in yard of 
Northeast Forestry University, at an interval of each two days 
from April 30 to June 2, 2004. The samples were quickly frozen 
with liquid nitrogen, and kept at -70 °C before extracting the total 
RNA. 
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Total RNA isolation 

The sample from the last time was used as origin of driver 
cDNA of Betula platyphylla, while the mixture from the other 
time was for the tester dDNA. 

The total RNA of Betula platyphylla was isolated by 
SDS/phenol method (Frederick et al. 1995).The integrity of total 
RNA was detected by the electrophoresis on a 0.8% agrose/EB 
gel, and its quality and amount were detected with 
UV-spectrometer (Eppendorf). 

cDNA synthesis 

First-strand cDNA synthesis was finished by the SMART™ 
RACE cDNA Amplification Kit. 

cDNA amplification 

PCR system includes 2.5 pL lOxLA buffer, 4 pL dNTP (2.5 
mM each), 0.5 pL primer (20 pM), 1 pL. First-strand cDNA, 1.5 
U LA-Taq polymerase (TAKARA), and sterile H 2 0 to total 25 
pL. 

PCR was performed at 94 °C for 1 min, then 94 °C for 30 s, 
68°C for 30 s, 72°C for 5 min for 18 circles, and with post incu¬ 
bation at 72’C for 5 min. The PCR product was purified with gel 
purification kit (Promega), and the 3 pL purified PCR product 
was detected by the electrophoresis of 1.0% agrose/EB gel. 

cDNA digestion by Rsa I, adaptor ligation to tester cDNA, 
and hybridization 

The procedure of Clontech PCR-Select™ cDNA Subtraction 
Kit was referenced. 

Amplification 

According to the procedure of Clontech PCR-Select™ cDNA 
Subtraction Kit, the diluted hybridization product was amplified 
for 25 circles with primer 1 by primary PCR, and diluted primary 
PCR product was amplified for 15 circles with nest primers in 
the second amplification. For T-A cloning, 15-min post incuba¬ 
tion was adopted at 72°C. 

After electrophoresis, small fragments shorter than 250 bp 
were discarded by gel cutting and the PCR product was purified 
with gel purification kit (Promega). 
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A subtracted cDNA library was constructed by the pGEM-T 
Vector System (Promega). After examination of ligation effi¬ 
ciency with M13 reverse primer and MI3 forward primer, a total 
of 150 randomly selected clones were sequenced to check for 
redundancy, which was 14%. 

Result 

Isolation of total RNA 

The integrity of RNA was examined by electrophoresing sam¬ 
ples on an agarose/EtBr gel. Total RNA typically exhibits two 
bright bands which correspond to ribosomal 28S and 18S RNA, 
respectively, with a ratio of intensities of about 1.5-2.5:1 (Fig. 1). 
As starting material, the quality of the total RNA can meet the 
needs for constructing a subtracted cDNA library by SMART 
strategy. 

Amplification of cDNA 

cDNA amplified with 18 circles was detected on an aga¬ 
rose/EtBr gel, and both driver cDNA and tester cDNA appeared 
smearing with a few bright bands (Fig. 2) which may be the 
transcripts of housekeeping genes. 



Examination of amplification of hybridization product 

Examined on an agrose/F.B gel, PCR product of fragments in 
the library appeared smearing with a size about 150-750 bp (Fig. 
3) in which the fragments were easy to be cloned to T-vector. 



Fig. 3 Amplification of hybridization product 

Lane M: DL-2000 marker; Lane I; PCR product 


Examination of ligation efficiency 

The average length of inserting fragments can be calculated 


according to the result from Fig. 4, which is about 350 bp, longer 
than the average length of Rsa 1 digestion (250 bp) because of 
gel excising. 



Fig. 4 Amplification of plasmids from random clones 

Lane M: DL-2000 marker 

Lane 1.3. 4: Products of recombinants 

Lane 2: 238 bp product of self-ligation. 


Analysis of the ESTs 

Among 150 clones, the ESTs (Expression sequence tags) with 
significant homology (E-value<le‘ 10 ) and sequences in the public 
databases were mainly classified into 5 clusters listed as follows. 


Cluster A: 

Suggested function 
Ribosomal protein 
PSII protein 

Ribulose-l,5-bisphosphate carboxylase large subunit 

ATPase proteolipid subunit 

ATP citrate lyase 

Isocitrate dehydrogenase 

Mitochondrial malate dehydrogenase 

Ubiquinone oxidoreductase 

Alcohol dehydrogenase 

Myo-inositol 1 -phosphate synthase 

Phosphoglycerate mutase 

Acetolactate synthase 

Acetyl-CoA carboxylase 

Gdsl-motif lipase 

Family II lipase EXL3 

Lipoxygenase 

Sphingolipid delta-8 desaturase 
RNA polymerase beta I subunit 
Glycyl tRNA synthetase 
ATP-dependent Clp protease 
Amine oxidase 1 
SSR alpha subunit 
Ubiquitin carrier protein E2 

26s proteasome regulatory particle non- ATPase subunit8 
Urease accessory protein G 


EST number 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 


The ESTs in cluster A reflect a background of gene expression 
in the process of floral development, which involves the process 
of photosynthesis, respiration, carbohydrates metabolism, lipid 
metabolism, and protein metabolism. Ubiquitin-dependent deg¬ 
radation of regulatory proteins can control many cellular proc¬ 
esses, including cell cycle progression, morphogenesis, and sig¬ 
nal transduction. Post-translational modification by ubiquitin 
targets many proteins for rapid degradation by the 26S protea¬ 
some. The E2 ubiquitin-conjugating enzyme is a component 
required for covalent attachment of ubiquitin onto lysine residues 
of the substrate (Galan & Peter 1999). Urease is a 
nickel-dependent metalloenzyme, and catalyzes the hydrolysis of 
urea to form ammonia and carbon dioxide. Urease plays an im- 
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portant role in germination and nitrogen metabolism of seedling. 
It may function coordinately with arginase in the utilization of 
seed protein reserves during germination (Sirko & Brodzik 2000). 
Urease accessory protein G is involved into urease maturation i.e. 
the placing of essential nickel in the urease active site. 


Cluster B: 

Suggested function EST number 

ABC transporter 26 

Nramp2 27 

Protein transport factor 28 

Amino acid transport protein 29 

Lipid transfer protein precursor 30 

Seven transmembrane protein Mlo4 31 

Membrane protein 32 

L-ascorbate oxidase 33 

Protein phosphatase 34 

Protein kinase 35 

BRI1-KD interacting protein 36 

NPG1 (No Pollen Germination 1) 37 

EICBP4 (ethylene-induced calmodulin-binding protein 4) 38 


Many ESTs in cluster B are related to the transmembrane 
transportation including ion, amino acid, protein and lipid. Three 
ATP binding cassette (ABC) transporter-like activities directed 
toward large amphipathic organic anions have recently been 
identified on the vacuolar membrane of plant cells. These are the 
Mg-ATP-energized, vanadate-inhibitable vacu olar accumulation 
of glutathione 5-conjugates (GS conjugates), chlorophyll catabo- 
lites, and bile acids, respectively. AtMRPl, from Arabidopsis 
thaliana that encodes a multi-specific ABC transporter compe¬ 
tent in the transport of both GS conjugates and chlorophyll cata- 
bolites (Lu et al. 1998; Liu et al. 2001). While protein kinase, 
BRI1-KD interacting protein, NPG1 and EICBP4 maybe is in¬ 
volved into signal transduction, especially NPG1 is a 
CaM-binding protein that mediates pollen germination and pol¬ 
len tube growth. It is expressed with a pollen-specific pattern in 
Arabidopsis thaliana , and is likely to be a key protein in the 
calcium signaling pathway in mediating calcium changes during 
pollen germination and tip growth (Golovkin & Reddy 2003). 


Cluster C: 

Suggested function EST number 

Cdc20 protein 39 

Microtubule associated protein 40 

Exostosin family protein-like 41 

Histone H 2 A protein 42 

Actin-related protein4 43 

Knolle 44 

SCC3 45 

Endo-l,3(4)-beta-glucanase 46 

Beta-glucosidase. 47 


The ESTs in cluster C are mainly related to cell cycle. His¬ 
tones are a major component of eukaryotic chromatin in associa¬ 
tion with DNA and nonhistone chromosomal proteins. These 
proteins and their mRNAs are usually synthesized in a cell cy¬ 
cle-dependent fashion at the beginning of S phase. This mode of 
synthesis is necessary to allow nucleosome formation on the 
duplicated DNA. Some histones are expressed in a non-cell cy¬ 
cle-dependent manner. In this case, mRNA of histone H 2 A is 
differentially expressed during development and that it exhibits a 
distinct cell-specific pattern of expression (Koning et al. 1991; 
Xu 1999). Actin is a ubiquitous component of the plant cy- 
toskeleton and participates in a number of important subcellular 
processes such as cell division plane localization, cell elongation. 


and cell shape determination. Actin is encoded by a gene family, 
and ACT4 is expressed in mature pollen and young vascular tis¬ 
sue in Arabidopsis (McKinney et al. 2002). SCC3 is a subunit of 
cohesions, protein complexes associated with sister chromatids 
in mitotic cell, and its expression should have something with 
regulation of chromosomal behavior (Pelttari et al. 2001). 
KNOLLE protein is related to a family of membrane anchored 
proteins, the syntaxins, which are implicated in directing intra¬ 
cellular vesicle trafficking. KNOLLE protein is specifically ex¬ 
pressed in mitotically dividing cells and in the cellularizing en¬ 
dosperm, mediating vesicle fusion in the plane of cell division, 
and was detected in mitotically dividing cells of various parts of 
the developing plant, including seedling root, inflorescence mer- 
istem, floral meristems and ovules, and the cellularizing en¬ 
dosperm (Lauber et al. 1997). Endo-1, 3 (4)-beta-D-glucosidase, 
beta-glucosidase is involved into the process of cell wall forma¬ 
tion and decomposition, is highly expressed during fruit ripen¬ 
ing. 


Cluster D: 

Suggested function EST number 

Manganese superoxide dismutase 48 

Metallothionein 49 

Cytochrome P450 50 

Polyphenol oxidase PPO1 51 


The ESTs in cluster D are mainly related to stress response. 
Expression of manganese superoxide dismutase appeared to be 
cross linked with the dehydration process during seed maturation. 
Metallothionein and cytochrome P450 can attach metal-ion, and 
play a role in decomposing toxin. PPO is an enzyme catalyzing 
the oxidation of o-diphenolic compounds to o-quinones as well 
as the hydroxylation of monophenols to o-diphenols. A role of 
foliar PPO in the defense against leaf-eating insects has been 
proposed and documented (Felton et al. 1989). Fruit commonly 
contains large amounts of PPO. However, the physiological 
function of PPO in fruit and other organs in healthy plants is still 
uncertain. 


Cluster E: 

Suggested function EST number 

DOF2 (DNA binding with one ftnger2) protein 52 

Nucleoid DNA-binding-like protein 53 

TFIID component TAF2 54 

Elongation factor 1 -alpha 1 55 

Eukaryotic initiation factor 4B 56 

GHMYB25 57 

SBP (squamosa promoter binding protein) 58 

HEN4 (Hua enhancer 4) 59 


The ESTs in duster E belong to regulatory factors including 
universal regulators and specific regulators. Dof domain proteins 
function as a transcriptional activator or a repressor were in¬ 
volved in diverse plant-specific biological processes and played 
critical roles as transcriptional regulators in plant growth and 
development (Yanagisawa 2002). TFIID component TAF2, 
elongation factor 1-alpha 1, eukaryotic initiation factor 4B and 
nucleoid DNA-binding-like protein are universe regulators, only 
nucleoid DNA-binding-like protein is localized in chloroplast. 
GHMYB25 is an MYB type transcription factor in Gossypium 
hirsutum that controls pigment formation, and there is a gene 
down-stream encoding isoflavone reductase in this subtracted 
cDNA library. SBP is a transcription factor regulating expression 
of API and Sqtia, floral organ identity gene in Arabidopsis 
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thaliana and Antirrhinum majus (Cardon 1997). HEN4 is a RNA 
binding protein in Arabidopsis thaliana that involved into RNA 
processing of floral organ identity gene AG (Cheng 2003). 

Discussion 

The quality of the subtracted cDNA library 

When total RNA acts as the starting material, a key factor af¬ 
fecting the quality of the subtracted cDNA library is the circle 
sum adopted in the process of cDNA amplification, which 
should meet the need of cDNA amount in subsequent procedure 
and avoid background as possible especially that from genes 
with high abundance. According to Clotech protocol, 17-19 is a 
reasonable sum range with starting material of 0.25-1 pg total 
RNA, and 18 is the optimized result in this report. 

Among the 150 clones sequenced, besides 25 fragments ho¬ 
mologous with sequences of unknown function in the public 
database, 8 shows no significant homology and 35 shows low 
homology, which seem to be due to short length, or due to large 
part of UTR. And another possibility is due to the genome iden¬ 
tity (some genes only exist in tree genome). Some ESTs in this 
cDNA library belong to regulators usually with low abundance 
in cDNA, which implicate a successfully uniformed process. 

In this study, gel excising was adopted to modify the result of 
blastX which might lost some information in the library, but it 
could be efficiency to screen valuable ESTs. 

Relationship between gene expression and development 

To enrich genes with differential expression in each develop¬ 
ment stage, this subtracted cDNA library is constructed accord¬ 
ing to temporal interval. To some extent, the sequencing result 
showed the background of gene expression throughout the period 
of sampling. During the 34 d of sampling, the female inflores¬ 
cence went through several development stages such as forma¬ 
tion of floral meristem, specification of floral organ, impregna¬ 
tion, seed maturation and final senescence. Expression of some 
genes {SBP, HEN4, NPG1 etc.) confirmed these events. 

Both API and AG are MADS box protein that play a dual role 
in floral development in Arabidopsis thaliana , i.e. API activates 
floral homeotic gene expression and specifies sepal and petal 
identity (Ng 2001), while AG controls floral meristem determi- 
nacy (Mizukami 1997) and specifies stamen and carpel identity. 
Expression of homology of SBP and HEN4 should be earlier 
than BpMADS3 and BpMADS6 (homolgy of API and AG in 
Betula pendula). It is indicated that their expression would be 
durable like BpSPLl did in Betula pendula reported before and 
AG did in Arabidopsis thaliana. It is required for analysis of 
function about the two regulators. 
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